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Abstract 

We present a pedagogical exposition of some applications of functional meth¬ 
ods in quantum field theory: we use heat-kernel and zeta-function techniques 
to study the Casimir effect, the pair production in strong electric fields, quan¬ 
tum helds at finite temperature and beta-functions for a self-interacting scalar 
field, QED and pure Yang-Mills theories. The more recent application to the 
UV/IR mixing phenomenon in noncommutative theories is also discussed in 
this framework. 


1 Introduction 

These lecture notes are intended to illustrate in a simple manner how spectral func¬ 
tions, viz. heat-trace and ^-function, are used to compute leading quantum cor¬ 
rections to physical quantities in field theory. The applications we consider are 
well-known and have been studied with one approach or another in most textbooks 
in Quantum Field Theory. Our main purpose is to readily provide, through a few 
plain examples, the basic tools of functional methods in held theory. 

Quantum effects in a particular held theory can be read from the spectrum of a 
certain differential operator. The effective action, for example, is given by the de¬ 
terminant on the inhnite-dimensional Hilbert space of an unbounded operator. This 
“inhnite” determinant can be written in terms of a divergent series. Fortunately, 
several criteria exist which allow us to unambiguously determine the sum of some 
divergent series. This procedure of assigning a hnite value to a divergent series is 
equivalent to the “regularization” of the UV-divergencies that occur in the computa¬ 
tion of one-loop Feynman diagrams. The functional method described in the present 
notes dehnes the value of a divergent series in terms of analytic extensions in the 
complex plane, the so-called “(^-function regularization”. The heat-kernel (or the 
corresponding heat-trace) is a mathematical tool closely related to the Qfunction 
which serves as a complementary regularization scheme. 

After introducing the heat-trace, the Qfunction and the main relation between 
them (section [2|) , we give an example where a (^-function is used to determine the 
value of some particular divergent series (section [S]) . To show a physical consequence 
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of this definition, we study the Casimir effect in its simplest version (sections 0] and 
ED- We also give an interpretation of this definition of divergent sums by comparing 
the (^-function method with the regularization in terms of a general cutoff function 
(section [6]) . Some relation between the use of analytic extensions and classical defi¬ 
nitions of divergent series (Cesaro, Abel and Borel summability) is illustrated with 
examples in appendix O 

Next we introduce one of the most powerful results in the theory of spectral 
functions: the asymptotic expansion of the heat-trace (section [7|); this expansion 
allows us to write a general expression for the effective action in terms of the Seeley- 
de Witt coefficients. As an example, we compute the /3-function of a real scalar 
field with a quartic self-interaction (section [H]) . The equivalence between the heat- 
trace and the (^-function regularizations is then understood from the point of view 
of renormalization theory in QFT (section [9|). 

The heat-trace can be used, in particular, to study the effect of quantum oscilla¬ 
tions of an electron/positron system in constant magnetic or electric backgrounds. 
In this context, we compute the /3-function of QED and the rate of pair production 
in the presence of strong electric fields (section [T0|) . 

Our last examples consist in the computation of the /3-function of Yang-Mills 
theories through the straightforward application of the heat-trace asymptotic expan¬ 
sion (section 1111) , of thermodynamical quantities in field theory at finite temperature 
(section fT^ and the propagator in a noncommutative field theory, whose nonplanar 
contribution exhibits the UV/IR mixing effect (section fTH|l . 

For the sake of clarity, we have chosen the simplest cases leaving aside any com¬ 
plexity which could obscure the presentation. Furthermore, the main mathematical 
results are stated without a detailed description of the assumed hypotheses, or a 
rigourous mathematical proof. As a compensation for these omissions we expect 
to provide -through quite explicit examples- the tools for a succinct calculation of 
these well-known quantum effects in field theory. 

2 Heat-trace and (^-function 

Heat-kernel and ^-function techniques deal with the calculation of the following 
related sums; 

and ^ A-% (2.1) 

n n 

where r > 0, s G C with 'R-{s) large enough and \n are eigenvalues of some positive 
definite (differential) operator, i.e., 

A4>n = Xn (pn with An > 0 . (2.2) 

In general, n represents a multi-index or even a continuous variable that characterizes 
the eigenstate. Of particular interest are the expansion of the first sum (heat-trace) 
for small r and the analytic extension of the second sum (('-function) to the whole 
complex s-plane. 

Let us begin by considering a second order differential operator A acting on a 
Hilbert space of functions "H; take, for example, the one-dimensional Laplacian 

A := -dl (2.3) 
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acting on the space % = ^^ 2 ( 0 ,1) of square integrable functions (/>(x) on the interval 
[0,1] with Dirichlet boundary conditions (/>(0) = (j){l) = 0. 

Let us then assume that T-L has an orthonormal basis of eigenfunctions 
(j)2{x), ... ol A with eigenvalues Ai, A 2 ,... In our example, the normalized eigenfunc¬ 
tions of the one-dimensional Laplacian 

4>n{x) = V 2 sin (nvrx) with n = 1, 2,3,... , (2.4) 

have eigenvalues 

An = , (2.5) 

which, as expected, grow to infinity with n. 

The heat-trace can now be defined as 

Tre-^^ = , (2.6) 

n 

for r > 0. In our example 

CX) 

^ , (2.7) 

n=l 

which can be written in terms of a Jacobi theta-function, 

CX) 

= ymmT)-\. (2.8) 

n=l 

Since the operator of our example does not have a zero mode, the heat-trace de¬ 
creases exponentially as r —)• -|-oo. More important to us is that, due to the infinite 
dimensionality of the Hilbert space, the heat-trace diverges for small r; the general 
behaviour of the heat-trace for small r is established in section (|7|) . 

The (^-function is defined as 

C(s) = TrAl-* = J^A-^ (2.9) 

n 

if TZ{s) is large enough, or as its (unique) analytic extension if not. In our example, 

C(s) = 7r-2^ Cr(2s) , (2.10) 

where C,r is known as the Riemann (^-function, and defined as 

00 

Ch(s) = E"’'- P.ii) 

n=l 

for lZ{s) > 1, or as its analytic extension otherwise. As a function on the whole 
complex plane Qr{s) has a unique simple pole at s = 1 and vanishes at negative even 
integers; B. Riemann conjecturecfl that any other zero of Cr{s) must be located at 
the line TZ{s) = 

^ “Hiervon ware allerdings ein strenger Beweis zu wiinschen; ich habe indess die Aufsuchung 
desselben nach einigen fliichtigen vergeblichen Versuchen vorlaufig bei Seite gelassen, da er fur den 
nachsten Zweck meiner Untersuchung entbehrlich schien” [T]. 
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We can now derive a very useful relation between both spectral functions. Ap¬ 
plying the relation (Mellin transform) 

1 r°° 

to each eigenvalue \n and performing the sum over the whole spectrum of A we 
obtain 

C{s) = r— dr Tr . (2.13) 

r(s) Jo 

This expression allows us to compute the analytic extension of the (^-function from 
the asymptotic expansion of the heat-trace for small values of r. 

3 Example: the Riemann (^-function 

As an example of the concepts of the previous section, let us consider the case in 
which the eigenvalues of the operator A are given by Xn = n = 1,2,3,... The 
corresponding (^-function is then the Riemann ^-function Cr{s), defined as the an¬ 
alytic extension of the series given by (|2.11D . which is convergent for TZ{s) > 1. In 
order to study its analytic extension to the whole complex plane we compute the 
corresponding heat-trace 


dr r® ^ e ^ 


( 2 . 12 ) 



n=l 


and use Mellin transform, given by (|2.13l) . to write 


Cij('S) 


1 


roo 

/ dr r' 

Jo 



(3.1) 


(3.2) 


As expected, the integral is well-defined only for TZ{s) > 1; otherwise it diverges at 
r = 0. In fact, the analytic extension of Cij('S) to the rest of the complex plane is 
determined by the behaviour of the heat-trace for small values of r. Therefore, we 
separate the integral at large values of r (say r > 1) from small values of r (then 
0 < r < 1) and use for the latter the expansion [2] 


r 


e’' - 1 


r- 1 lr2 

^ " n! 2^^ 6 2! 

n=0 


^ 42 ¥ 


(3.3) 


which is actually valid for |r| < 27r; Bn are the Bernoulli numbers. After integrating 
this powers series in the interval [0,1] we obtain 


Cr{s) 


r(s) ^e^-l + 


(3.4) 


Although the previous manipulations are only valid for Tl{s) > 1, expression (|3.4p 
now gives the analytic extension of Cr{s) to the whole complex plane: the first 
term between brackets is an entire function whereas the series in the second term is 
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analytic for all s G C, except where it presents simple poles -namely, at s = 1 and at 
all negative integers. However, since the term l/r(s) vanishes at negative integers 
then the Riemann ^-function has a nnique pole at s = 1, the residue being Bq = 1. 

The term l/r(s) also facilitates the determination of Cr{s) at each negative 
integer, which is given by the corresponding residue of the expression in brackets. 
For instance, at s = 0 only the term in the series corresponding to n = 1 gives a 
nonvanishing value, 

CK(0)=Hi = -i. (3.5) 

Since the value of Cr{s) at a negative integer —k is proportional to and 

Bernoulli numbers of odd order vanish -except from Bi- then the Riemann (^- 
function vanishes at even negative integers; these are the trivial zeroes of Cr{s)- 
On the contrary, at odd negative integers + 1) = —B2k/‘2k, so that, e.g., 

&(-!) == (3.6) 

In summary, expression (12.lip is only valid for TZ{s) > 1 but it suffices to define 
Cr{s) in the whole complex plane as its unique analytic extension. In the (^-function 
regularization scheme, the analytic extensions we have just computed are used to 
define the sum of the following divergent series: 


1 + 1 + 1 + .. 

• := C/?(0) = , 

(3.7) 

1 + 2 + 3 + .. 

• := Ck(-1) = , 

(3.8) 

1 + 4 + 9 + .. 

. :=Cij(-2) =0. 

(3.9) 


Certainly, any definition might seem uncanny insofar as it assigns a vanishing or even 
a negative value to a sum of positive terms. In section [6] -where we regularize these 
divergent series with a smooth function of compact support- we will understand 
how these unexpected values appear after subtracting the infinite contribution orig¬ 
inated by the removal of the cutoff. In appendix!^ a comparison between analytic 
extensions and other definitions of divergent sums is sketched. 

4 Application: vacuum energy of a particle in a ring 

Before giving an insight into this intriguing regularization of divergent sums, we 
apply the results they provide to the calculation of a physical quantity, namely the 
vacuum energy of a ring generated by the quantum oscillations of the field of a 
scalar particle. Later, in section [5l we will compute the vacuum energy in a more 
realistic setting where the predictions of the ^-function regularization have been 
experimentally confirmed. 

Let us consider a massless scalar particle on a ring of radius R. Its dynamics in 
the context of relativistic quantum mechanics is described by a wave function (f){x,t) 
that satisfies 


{dt - t) = 0 . 


(4.1) 
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This equation has infinitely many stationary solutions which can be identified with 
an integer re G Z, 




ipnX-iUnt 
^ 1 


where 


Wn 



(4.2) 


(4.3) 


and Pn = n/R, due to periodicity in 2'kR. From the point of view of classical 
mechanics, each of these solutions with momentum n/R and positive energy is 
a travelling wave that can also be considered as a normal mode of frequency ujn 
of a system of infinitely many coupled harmonic oscillators. The quantization of 
the field t) is equivalent to the quantization of this infinite set of harmonic 
oscillators. As is well-known, due to quantum oscillations, the minimum energy of 
a single harmonic oscillator of frequency lo is given by As a consequence, the 

minimum energy of the quantized field should correspond to the (infinite) sum of the 
ground state energies of all the normal modes of the system of coupled oscillators. 

In the language of quantum field theory, the space of physical states of the 
quantized field is spanned by vectors characterized by the number of scalar particles 
in each of the stationary solutions (|4.2I) . In particular, there is a vector |0) which 
represents the vacuum state, i.e., the state that describes the system in the absence 
of such particles. Surprisingly, although the total number of particles in the vacuum 
state (0|A^|0) vanishes, its energy Eq = (0|i?|0) is not necessarily zero. As explained 
above, in the case of the particle in a ring, Eq corresponds to the infinite sum 



(4.4) 


which -interestingly- depends on the size of the ring. As a matter of fact, it is obvious 
from dimensional reasons that Eq should be proportional to 1/R but we obtained 
for the proportionality coefficient a divergent series of positive terms; should it be a 
positive number? Is there a way to get some physical information from expression 
(14.4p ? In the ^-function regularization described in the previous section, the vacuum 
energy Eq is defined as 


^0 = ^ C/?(-l) 


1 1 
12 R ■ 


(4.5) 


This result indicates that, even in the absence of the massless particle, the quantum 
oscillations of its vacuum state originate an attractive force under which the ring 
tends to shrink to zero siz^. 


5 Casimir force between two conducting plates 

In this section we will determine the Casimir force between two parallel plates; our 
purpose is to give some more details of the basic framework of quantum field theory, 

^This force could be considered as responsible for the smallness of compactified extra dimensions 
in Kaluza-Klein models [3]. 
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as well as to illustrate some usual manipulations regarding the computation of the 
analytic extension of the i^-function in terms of the asymptotic expansion of the 
heat-trace. 

Let us then consider a relativistic massive particle confined between two plates 
separated a distance L in the x-direction. In the context of relativistic quantum 
mechanics the particle is described by a wave function (/)(x, y, t) -with support in 
0 < X < L, y £ and t G M- that minimizes the action 

'S’[0] ~ [ ^^ [ k . (5.1) 

«/ M J J 0 

These configurations therefore satisfy 

(52 + A)</. = 0, (5.2) 

where A is the positive, second order differential operator in 3-dimensional space 
given by 

A := —d^ + w? . (5.3) 


In addition, we impose some, say Dirichlet, boundary conditions at x = 0 and x = L 
that represent the effect of the plates. Let us compute the spectrum of this operator. 
The normalized eigenfunctions of A are given bjo 


y) = \Ij {n-Kx/L) , (5.4) 

where G and n G Z"*"; the corresponding eigenvalues can be written as uj‘i , 

k,n 

with 

= Y • (5-5) 

In consequence, stationary solutions to eq. (j5.2p can be written as 

(^-6) 

Now, we are ready to obtain from the spectrum of the operator A some physical 
information about the quantized field. 

In quantum held theory, the scalar held is given by an operator built from a 
superposition of the normal modes (15.61) . 




®These eigenfunctions satisfy the following orthogonality and completeness conditions: 

f <fy f dx 4>l^{x,y)(j}p ^,{x,y) = {27vf 5{k-k')5„„', 

Jo 

r k 

) _ 
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Now, equal-time canonical commutation relations between (f){x, y, t) and its conju¬ 
gate field lead to the following algebra of creation and annihilation operators: 

This algebra can be represented in the Fock space generated by the repeated action 
of creation operators at on the vacuum state |0), which is a vector annihilated by 

all H,n 

Moreover, from the classical action m one obtains the Hamiltonian of the 
quantized field, which corresponds -if written in terms of creation and annihilation 
operators- to an infinite set of uncoupled harmonic oscillators (the normal modes 
described in section S]). In fact, the vacuum expectation value of this Hamiltonian 
(per unit area in the directions parameterized by y € results 


Eo = {0\H\0) 


If £k ^ 

2 L (2vr)2 A- 


(5.9) 


which, as mentioned in section [H is the sum of the ground state energies of an 
inhnite set of harmonic oscillators with frequencies 

Expression (15.9p is of course ill-defined since the series and the integral are not 
convergent, due to the behaviour of for large n and \k\', this is a manifestation 
of the UV-divergencies one encounters when compnting qnantum corrections in field 
theories. The (^-fnnction method provides a regnlarization of this divergence by 
means of the following function of the complex variable s G C: 


Eo{s) 


y f d?k 
2 ^2^ 


E 

n=l 


U 


k,n 

k- 


-2s 


2 


C(s) 


(5.10) 


where C('S) is the ^-function of the operator A, defined as 

cPk 






.2 s 
k,n' 


(5.11) 


for 'R,{s) > |, or as its analytic extension otherwise. Note that, for dimensional 
reasons, we have also introduced a parameter /r as an arbitrary scale with mass di¬ 
mensions whose physical significance will be elucidated later. Finally, by comparison 
with eq. (15.9p . the vacuum energy is defined as the analytic extension of E'o(s) to 
s = — ^. If the (^-function were analytic at this point of the complex plane then the 
vacuum energy would be -as one would expect- independent of the arbitrary scale 
y, but this might not be the case. 

As we did in section [3] for the Riemann ^-function, in order to obtain the analytic 
extension of (|5.11l) we first compute the heat-trace of the operator A, 


Tie-^^ 


f (fk 

L'iEf 


E 

n=l 


L 


(fk 

R2 (2vr)2 


E 

n=l 



e-Tvn? 

(dvrr) 2 


L — -|- 2L 



(5.12) 










where, after integrating in A: G we have rewritten the sum by using Poisson 
inversion formula, 


E»- 


E' 

n^’L 


(for a > 0), 


(5.13) 


because it makes explicit the divergent behaviour of the heat-trace for small values 
of the parameter r. Prom relation (j2.13h we obtain for the (^-function 


C(s) = 


1 


r(5: 


f 


dr r 


S —1 


(dvrr) 


L — yE r 2 -|- 2L ( 


71=1 


Pj-g _ 3^ 

(47r)2 r(s) 


m? 1 m2 2 .^2^ K3_^{2mLn) 

Stt s —1 27rtr(s) n=i 


(5.14) 


It is important to remark that the terms in the heat-trace that decrease exponentially 
as r —)• 0"*“ give contributions to the ('-function that are analytic in the whole complex 
plane. In conclusion, ((s) has a simple pole at s = 1 (due to the second term in 
(15.141) 1 and an infinite number of other simple poles at s = | — Z+ (due to the 
first term in (j5.14p l. In particular, it has a simple pole at s = — ^ with residue 
—m^L/327r^. The function Eq{s) can thus be written at s = — ^ e as 


mlL 1 
647r^ e 


Eo{-\ + e) = fin(m,/r, L) 0(e), 


where the finite part reads 

m^L 




fin(m,/r,L) = 77 :^ + 


247r 6477^ 


/ m?\ 1 

log ^ -h - - 2 log 2 


\E 


m 

Sit'll 


E 

n=l 


(5.15) 

K2{2mLn) 




(5.16) 


As was already seen in expression (I5.1UI) . a pole in ((s) at s = — ^ brings forth a 
(logarithmic) dependence of the vacuum energy on the arbitrary scale fi. In the 
massless case, where the only dimensionful parameter is L, the analytic extension 
of Eq{s) is finite at s = — ^ and thus /r-independent; therefore, the limit e —>■ 0 
of expression (I5.15P gives an unambiguous result for the vacuum energy (per unit 
area), 


= (form = 0). (5.17) 

However, the simple pole of Eo{s) at s = — ^ for the massive case requires some 
additional physical prescription to define a renormalized vacuum energy. For in¬ 
stance, if the field is infinitely massive then quantum oscillations are blurred and, 
consequently, the vacuum energy is expected to vanish; thus an admissible procedure 
consists in removing (by a minimal subtraction) those terms in eq. (I5.15p which do 
not vanish in the m —)■ 00 limit. After this subtraction we obtain for the vacuum 
energy of the massive field 


En — — 


m 

Sit'll 


E 

n=l 


K2{2mLn) 




(5.18) 
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Figure 1: Casimir attractive pressure (in mPa) as a 
function of the plates separation (in /rm) for mass¬ 
less and massive (m = 1 eV) fields. 

Taking into account that the physical consequences of the vacuum energy are associ¬ 
ated with its variation with the separation of the plates, one could also remove from 
expression (I5.15P those terms which are independent of L, as well as those linear in 
L -which together with the vacuum energy in the space outside the plates give an 
T-independent contribution of the vacuum energy of the whole space The result 
obtained with this prescription is again given by expression (I5.18p . 

In conclusion, there must exist an attractive force (per unit area) between two 
Dirichlet-type parallel plates, that is originated by the quantum oscillations of the 
vacuum state of the field. Since the force is exponentially decreasing for massive 
fields, the most prominent contribution corresponds to the oscillations of massless 
particles. The pressure on the parallel plates is given by the variation of the vacuum 
energy with their distance, which for the massless case reads (see eq. (j5. 171) 1. 

p-.= -dLEo = -nc^^. (5.19) 

This pressure is represented in figure [H as well as for the massive case (see eq. 
(|5.18p l. as a function of the separation of the plates. As the figure shows, at a plates 
separation of 1 /iin, massless fields originate a small pressure of the order of mPa; for 
a mass of 1 eV the pressure decreases to 0.3% this valutEI, whereas it is completely 
negligible for masses of the order of MeV. 

This effect was predicted in 1948 by H. B. G. Casimir [4], who studied the attrac¬ 
tion experienced by two parallel conducting plates due to the vacuum oscillations of 
the electromagnetic fielq^. Already in 1958 this Casimir force could be experimen¬ 
tally observed [5]. Since then, much more accurate experiments have successfully 
measured the Casimir force in different settings PE] and are used nowadays to test 
many theoretical models that take into account more realistic aspects, such as the 
finite conductivity of the boundaries, the roughness of their surfaces as well as the 
effect of temperature. 

"^Nevertheless, this exponential decrease of the Casimir force with the distance for massive fields 
is not verified, in general, in the presence of curved boundaries. 

®The vacuum pressure originally derived in [4] is twice expression (ITT9t because of the two 
transverse polarizations of the electromagnetic field. 
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6 Cutoff regularization 


In the previous sections we have dehned the result of a divergent series in terms of 
the analytic extension of the ^-function and we have found that this result leads to 
physical predictions -the Casimir force- that have been conhrmed by very accurate 
experiments. Nevertheless, a physicist might be more familiar with the regulariza¬ 
tion of a divergent sum in terms of a cutoff functiorj^. In this section we will compare 
both procedures in order to become more familiar with the meaning of the dehnition 
in terms of analytic extensions. 

Based on Euler-Maclaurin formula 
N 


iV „7v 

~ / * 
n=l ■^0 


dxfix) + I IfiN) - /(O)] + E /(2-i)(iV) - /(2-i)(0) 




- (2n)! 1' 


( 6 . 1 ) 


Ramanujan computed the value of a divergent series using an appropriate cutoff 
function [8]. We will depict his procedure by computing the same inhnite sums 
which we have already computed using the (Riemann) ^-function regularization. 
Let then uj : —)• M be a smooth, positive and bounded function with support in 

the interval [0,1] and such that, for some e > 0, takes the constant value lj{x) = 1 
for any 0 < a: < 1 — e. If we apply Euler-Maclaurin formula to the function f{x) = 
x^ uj{x/N), for k G Z+, we readily obtain 


poo TD 

Y^n^uj{n/N) = / dx x^oj{x) ^ , ( 6 . 2 ) 

n=l -^0 ^ ^ 


where the lower sign should be used only for the particular case /c = 0. We conclude 
that if we replace the infinite sum of the A;-th powers of positive integers with the sum 
of the hrst N terms -smoothed with a cutoff uj{n/N)- we obtain an expression which, 
unsurprisingly, diverges as N —)• oo. However, this divergence can be isolated in a 
term which corresponds to the integral in Euler-Maclaurin formula (proportional to 
the fe-th moment Ck[oj\ of uj{x)). The remainder is a hnite term which is independent 
of the cutoff and coincides with the value given by the analytic extension of the 
corresponding (^-function; Ramanujan claimed this term “is like the center of gravity 
of a body”. In particular, for fe = 0,1, 2, we obtain 

OO ^ 

^uj{n/N) cq[uj\N - (6.3) 

n=l 

oo ^ 

^ nu)(n/N) ~ Cl [w] > (6-4) 

n=l 

oo 

^ w(n/A^) ~ C2[w] . (6.5) 

n=l 


Note that after subtracting the leading (divergent) behaviour one obtains the values 
of the Riemann (^-function Cr{s) at s = 0,1, 2. Erom this point of view, the analytic 
extension Ci?(s) gives the difference between the divergent “smoothed” series and 
the corresponding “smoothed” integral. 

® Actually, H. B. G. Casimir used a cutoff function in [4], whose physical meaning he interpreted 
as the fact that for high energy virtual photons the conducting plates must become transparent. 
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7 Asymptotic expansion of the heat-trace 

Now we will state a very important mathematical result that describes the asymp¬ 
totic expansion of the heat-trace Tre“’^"^ for small values of r. We have already 
seen that this expansion allows the computation of the analytic extension of the 
^-function; in particular, it provides complete information of its pole structure. 

Let us consider a second order matricial differential operator A acting on the 
space of n-component functions C”(8)L2(1^^)- Such an operator defines a generalized 
heat-equation (or evolution equation in Euclidean time, if A were a Hamiltonian) 

(dr + A) x) = 0 (7.1) 

whose solution -under some given initial condition i^(0; x)- can be written in terms of 
the integral operator called heat-operator (or evolution operator in Euclidean 

time), as follows, 

(j)(T]x) = e~'^^4>(0]x)= f dx K(t-,x,x') (j)(0;x') . (7.2) 

Jm.d. 

If, as before, C" (8) L 2 (M'^) has an orthonormal basis of eigenfunctions ■ ■ ■ of 

A with eigenvalues Ai, A 2 ,..., then the heat-kernel K(t\ x, x') can be written as the 
series 


CX> 

K{r; X, x') = ^ (j)n(x) (t)l,(x '), (7.3) 

n=\ 

and the heat-trace is 

„ 00 

Tre“'^"^= / dxtr K(t;x,x) = 'y^ , (7.4) 

where tr represents the finite trace in the C” part of C” (8> L 2 (M'^). 

As the proper time r tends to zero the heat-operator tends to the identity op¬ 
erator in an infinite dimensional Hilbert space, so the heat-trace diverges in this 
limit. The asymptotic behaviour of the heat-trace for small values of the proper 
time has important applications in QFT and has been systematically described for 
a quite general class of differential operators defined on smooth manifolds. It has 
been proved that, under quite general conditions, the following asymptotic expansion 
holds [9]: 


Tr e 


—r A 


(dyrr)'^/^ _ 


On (A)- 


n=0 


(7.5) 


where the Seeley-de Witt coefficients an(A) are integrals of local expressions of the 
coefficients of the differential operator A (and the geometric properties of the base 
manifold in the case of curved spacetimes). If the base manifold has boundaries then 
the sum in (j7.5l) also includes positive half-integers n G Z'*'/2. Nevertheless, in our 
applications we only need to consider the simpler case of the flat spacetime without 
boundaries 
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Let A be the second order differential operator 


A = -(lndi+CJi(x)f + V(x), 


(7.6) 


where uJi{x),V{x) are matrix-valued functions that take values in C"' (8> C”. Then 
the first Seeley-de Witt coefficients of expansion (17.51) are given by the following 
expressions: 


ao{A) = 

J dx tr In , 

(7.7) 

ai{A) = 

— / dx tiV(x), 

(7.8) 


a 2 {A) = / dx tr {^V‘^{x) - I {di+u;ifV{x) + - dju;i + [uji,u!j]f} . 

(7.9) 


Since we are considering a non-compact base manifold, these integrals are in general 
ill-defined; one should consider instead the expressions corresponding to (|7.5p -as 
well as to (1131), dZSl) and dZS])- for the local quantity K{t;x,x) or, alternatively, 
for the trace of the heat-operator but including a convenient smearing function. Yet, 
as long as we appropriately take into account the infinite volume contributions, we 
will be able to derive from these expressions the effective action of a field theory 
over the whole 

The asymptotic expansion of the heat-trace, together with relation (I2.13p . allows 
us to determine the pole structure of the (^-function by the same procedure followed 
in section [3] for the particular case of the Riemann ((^-function, 


E 

^ n=0 



1 

(dyrr)'^/^ 

an(A) 

I 

(47r)<i/2 

s — {d/2 


n=0 


-|- (analytic terms) 


(7.10) 


We conclude that, in general, the ^-function has simple poles at d/2 — -except 
at negative integers, due to the cancellation with l/r(s)- and that the residues are 
given by the Seeley-de Witt coefficients. Moreover, the value of C{s) at a negative 
integer s = —k € Z“ is given by 


a-k) 




«d/2+fc(^) 
(47r)<^/2 ’ 


(7.11) 


for even d, but vanishes in odd dimensions. 


8 The effective action 

In a field theory described by the classical action 5 [(/>], the effective action r[(/)] 
is another functional of the field cf) which -according to Feynman’s quantization 
procedure- is obtained after averaging all field configurations ip with the weight 
g-S[i,p]_ gince S'):/?] is measured in units of h, the leading contribution to the effec¬ 
tive action in an expansion in powers of h is given by the classical action. In this 
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section we will show that the first quantum correction is given by the functional de¬ 
terminant of a differential operator determined by the second functional derivative 
of the classical action. We will then describe how the spectral functions are used to 
compute this functional determinant and we will apply the asymptotic expansion of 
the heat-trace to characterize the UV-divergent terms in the effective action. With 
this procedure we will finally compute the /3-function of the self-coupling constant 
of a scalar field. 

We begin by defining the partition function Z[J] and the generating functional 
W[J] -both depending on some arbitrary source J{x)- as 

Z[J] := [ Vip 3 --^ . (8.1) 


The mean held (also understood as (0|(^(x)|0)) in the presence of this source is 

the functions J{x) and </>(x) are thus called conjugate helds. Finally, we dehne the 
effective action as 

m-=W[J] +I J<p, (8.3) 


where J must be understood as implicitly determined by cf) through the relation 
()8.2p . From this dehnition it immediately follows the equation 

= Jix) 

In consequence, in the absence of the “external source” J{x), the vacuum expec¬ 
tation value (j){x) of the held minimizes the effective action F, unlike the classical 
conhguration of the held, which minimizes the classical action S. Moreover, if we 
make a functional expansion of the effective action in powers of the held (j), the 
coefficient of the n-th power of the held provides the (proper) scattering amplitude 
of n particles. In this sense, the effective action describes the full behaviour of the 
system including all quantum effects. 

Now it is convenient to consider (j){x) as an arbitrary conhguration -determined 
by an appropriate J{x)— and to shift the integration variable as (p ^ (j) + (p io obtain 

Z[J] = y Py, , (8.5) 

where 5^S denotes the functional derivative 5S[ip\/5ip{x) evaluated at the conhgu¬ 
ration (/>(x). Similarly, the kernel 5‘^S is given by 5“^S[p\/5‘p{x)5‘p{x') evaluated at 
(l){x). The effective action then reads 



(8.4) 


r[(/)] = 5[(/>] - h\og j Vip e s^r)v^+Vno{v’3) ^ 


where we have rescaled (p —)• y/h(p. Since F = S + 0{K) we can write 

r[(/.] = s[^]-hiog J Vip (i + o(n)), 


(8.7) 
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where 0{h) comes from expressions which are quartic in (p. The Gaussian functional 
integral is given by the functional determinant of the operator A which defines the 
quadratic form in the exponent, that is, 

T[cj)] = S[^] +^logDet A+ (8.8) 

This expression gives the leading quantum contributions (one-loop corrections) to 
the effective action. The operator A, sometimes referred to as the operator of quan¬ 
tum fluctuations of the field, is defined by the kernel 6 ‘^S. Take, for example, the 
action of a real scalar field cj) on (Euclidean) with a quartic self-interaction 

S[(p]= f dx ^ ^ X . (8.9) 

Jr* 

Then the operator of quantum fluctuations is the second order differential operator 

A = —d^ + ^ X(j)‘^ . (8.10) 

If we assume that A has eigenvalues Ai, A 2 , A 3 ,..., which in general grow to infinity, 
we must now indicate how to compute its determinant. The (^-function and the 
heat-trace provide different definitions of the functional determinant. In principle, 
we would expect 


log Det A = ^logAn, (8.11) 

n 

but this series is in general divergent. In the (^-function approach one introduces 
an appropriately decreasing power of the eigenvalues A“^ and then computes the 
analytic extension to s = 0 of the regularized series. 


logDet^ := ^log(An//i^) (A„//r^) * . (8.12) 

n s=0 

The series is convergent for 'TZ{s) large enough and the determinant is defined as the 
analytic extension of this series to s = 0. As before, we have introduced an arbitrary 
parameter p with mass dimensions whose relevance in the determination of physical 
quantities will be discussed later. This definition can also be written as 

logDet A =-C'(0) - C(0) log/i^ , (8.13) 


where 


n 


(8.14) 


is the C-function of the operator A. As we have seen from (I7.10p . the ^-function is 
analytic at s = 0. 

The heat-trace instead provides a different definition, based on the identity 




r 



log(AJA2)-7 + 0(An/A2), 


(8.15) 
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where 7 is Euler-Mascheroni constant and A is some high-energy cutoff which could 
be eventually removed by making, when possible, A —)• 00 . In this approach the 
determinant is then dehned as 

roo j 

logDetA:=-/ —Tre"^^. (8.16) 

iA-2 r 

In this case, the arbitrary mass scale is introduced by A. The limit A —>■ 00 charac¬ 
terizes the UV-behaviour of the theory. 

Let us now compute the one-loop effective action for the self-interacting scalar 
held in described by the classical action ()8.9p . If we use the heat-trace dehnition 
of the functional determinant we obtain 

= S[^] - - r — e-™' Tr , (8.17) 

2 Ja-2 t 

where we have omitted h. Note that we have factorized the mass term in A (see 
() 8 . 101) 1 from the heat-trace. The function V(x) is the held-dependent expression 

V{x) = ^X^‘^. (8.18) 

Next, we insert in (18.171) the heat-trace expansion as given by expression (17.5p . 
together with the expressions (|7.7p . (17.8p and (17.9p for the coefficients, to obtain 

1 Ht 9 1 

r[(/>] = S[4>\ - - — e"™ ——^ {ao + aiT + a 2 T‘^ + ...} 

2 t (dvrr)^ 

=- 3 I? C ^ {I" I^ A" T ■ ■ ■ } ■ 

(8.19) 

In the last line we have only displayed the three terms that lead to divergencies in 
the effective action as A —)• 00 : a held-independent term, a mass term (quadratic 
in the field) and a quartic self-interaction term. We therefore face two different 
problems: the dependence of the effective action on the cutoff A and the inhnities 
that appear when the cutoff is removed. Actually, these two unacceptable results 
are related and can be circumvented on the basis of a same assumption. 

Note that -taking into account the cosmological constant- these three types of 
terms are already present in the classical action. The renormalization procedure 
consists in regarding the original action S in eq. (18.9p . as well as its parameters m, A 
(and the cosmological constant), as the result of integrating out (with some unknown, 
eventually more complex action) held conhgurations which, in Fourier space, have 
momenta above the arbitrary cutoff A. Thus, the action in expression (18.9p only 
describes held conhgurations with momenta below A. Prom this point of view, the 
parameters m, X clearly depend on A and it is to be expected that this dependence is 
such that the effective action -from which physical results are obtained- is hnite and 
A-independent. Accordingly, the effective action can be written in terms of physical 
parameters a^ 

r[(f)] = [ dx {l{d(f)f + I m^hys + i Vys > (8-20) 

^In this model, after considering two loop corrections also the field (j) must be renormalized. 
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where ^phys Aphys (of course, A-independent) are determined by physical mea¬ 
surements. Comparing expressions (I8.19P and (|8.20l) we readily obtain 


2 2 
"iphys = rn 




= m 


1 + 


A 

327r^ m? 



\phys 


= Aa- 


3A 

327r2 


dr 


' vn?-1 


= \{\- 


3A 

IGvr^ 


log {kim) + ... 


( 8 . 21 ) 

( 8 . 22 ) 


These expressions intrinsically determine the dependence, up to one-loop correc¬ 
tions, of the original parameters m, A on the cutoff A. The /3-function, given by the 
derivative of the coupling constant with respect to the cutoff A, then reads 

/3(A) :=A5aA = ^A2. (8.23) 

This expression is the leading quantum contribution; higher order quantum correc¬ 
tions contain higher powers of A. Therefore, as long as the coupling constant is small, 
since its derivative is positive, it becomes larger as A increases or, equivalently, as 
the theory involves field configurations with higher momenta. This justifies the use 
of perturbation theory to study the scattering amplitudes of these self-interacting 
scalar particles at low energies. However, if the coupling constant increases to an 
infinite value at some finite A then the full theory becomes inconsistent, unless it is 
trivial; this is called Landau pole problem, and is also present in QED (see section 

uni). 


9 Equivalence between both regularizations 

As we have seen, the (/"-function and the heat-trace give two different definitions of 
the functional determinant. Let us briefly study how their difference arise in the 
computation of the effective action for a field theory described by an operator of 
quantum fluctuations A. 

As shown in the previous section, it is convenient to separate the mass term from 
the operator A so we assume in this section that the second functional derivative 5‘^S 
is the kernel of the operator A -|- m?. The functional determinant in the heat-trace 
approach reads 

poo —rm? 

logDet [A + rm?) = — dr - Tr 

- CXD 

=-713372 «„{. 4 ), ( 9 . 1 ) 

^ ' 71=0 

in terms of the incomplete gamma function. Let us analyse the UV-divergent terms 
in this series separately. For 0 < n < d/2 we get the surcH 

d/2 

-33372 E r (-1 + ". 5l) “niA) , (9.2) 

' ' n=0 

^Strictly speaking, what follows holds for even d; for odd d one must consider 0 < n < d/2 — 1/2. 
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which diverges as for A —)• oo; the term corresponding to n = (i/2 (only 

present in even dimensions) gives a logarithmically divergent contribution. These 
divergencies are removed by introducing similar counterterms in the classical action: 
if some of these terms are already present in the action then the corresponding 
constants get renormalized; otherwise, one says that new couplings are generated by 
quantum corrections. In any case, the theory can be made one-loop finite if the action 
contains or admits terms of the form given by ao(^), ai(A), 02 ( 2 !),... 

Once the divergencies in these terms are removed by a redefinition of the corre¬ 
sponding parameters in the classical action, the remaining terms in the functional 
determinant are given by 

r ( n - I - 1 ) 

- ad/2+i+n{A) (for even d) (9.3) 

- ^ ad/ 2 +i/ 2 +n{A) (for odd d) . (9.4) 

' ' n=0 

Let us now study the determinant as defined in the ^-function approach, 

logDet (A + m^) = —C'(0) — C(0) log . (9.5) 


By means of the Mellin transform the ^-function can be written as 

1 

drr^ " e 

d- 2 n^{s - d/2 + n) 


1 r 

= fM / 

r(s) Jo 

— 2 <? ^ 
m 

(47r)'i/2 


(47rr)<^/2 
^ ' n=0 


m 


n=0 


r(s) 


an{A). 


(9.6) 


Let us consider first the value of (/(s) at s = 0. For odd d the term l/r(s) cancels 
any contribution. In consequence, in odd dimensions (/(O) = 0 and the determinant 
is independent of the arbitrary scale /r. In even dimensions instead, the term l/r(s) 
cancels the divergencies from the terms with 0<n<d/2ats = 0 and the only 
non-vanishing contributions are 


C(o) 


(- 1)^/2 ^ (-l)n^d- 2 n 

( 47 r )‘^/2 ^ {d/2-ny. 


an{A ). 


Let us now compute the derivative of the C-function, 


(9.7) 


- C'('S) = logm^C('S) - 

- TTW £ + ") - 'HO] , (9.8) 

(4vr)'^/2 ^ r(s) 

in terms of the gamma and digamma functions. At s = 0 the first term is propor¬ 
tional to C(0) and together with the second term in the r.h.s. of (19.5p gives a total 
contribution log {m?/fj?) ("(0) to the functional determinant. In consequence, the de¬ 
pendence on the arbitrary scale /r is proportional to the coefficients ao(A), ai(A),..., 
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ad/ 2 {^) aiid is removed from physical quantities by introducing an appropriate de¬ 
pendence on the scale // of the original parameters in the action. This is equivalent 
to the dependence of the parameters on a renormalization scale // when imple¬ 
menting renormalization prescriptions in terms of scattering amplitudes at some 
arbitrary scale of external momenta. On the other hand, the terms in C^(0) with 
0 <n<d /2 — 1 give contributions to the effective action of the same form as the 
ones given in (IHZD; thus they can also be considered as part of the counterterms 
that renormalize the classical actiorH. 

As regards the remaining terms in the series in (19.8p . in even dimensions the 
term n = d/2 vanishes, whereas the terms with n > d/2 give the series 

~ U^YI2 X] ®rf/2+l+n(^), (9.9) 

' ' n=0 

which coincides with the heat-kernel result for the case of even d (see (19.30 . Anal¬ 
ogously, on can immediately prove that for odd d we also get the same result as in 
the heat-kernel regularization. 

In conclusion, both in the heat-trace and in the ^-function regularizations one 
introduces arbitrary scales, A and /x, respectively. A finite number of terms of the 
effective action diverge with the cutoff A or depend on the arbitrary scale their 
dependence with the field is, in both cases, given by the first Seeley-de Witt coeffi¬ 
cients On, with 0 < n < (i/2. Therefore, the difference between both approaches lies 
in the dependence of the coupling constants with the large cutoff A or the arbitrary 
scale /i. This dependence is established through the experimental measurement of an 
appropriate number of scattering amplitudes; once this dependence is determined, 
the remaining contributions to the effective action (see eqs. (|9.3I1 and (19.41) ') are finite 
and unambiguously given by the Seeley-de Witt coefficients of higher order. 

10 Electrons in a constant magnetic field 

In this section we consider the coupling of electrons to a constant magnetic field 
and study, to leading order, the effect of the quantization of the electron field. Since 
fermions are quantized in terms of Grassmann variables, the one-loop contribution to 
the effective action of an electron/positron system in an electromagnetic background 
readJ^ 

r = S-log j e- ^ = S - logBet {0 - m} , (lO.l) 

where S is the classical action for the electromagnetic field. We will use the heat- 
trace to compute this functional determinant. However, this method applies for 
positive operators so we need to perform another formal manipulation. Since in four 
dimensions there exists a charge conjugation matrix C which satisfies C'y^C~^ = 
— ( 7 ^)”^, we can write 

logDet {il/) — mj = logDet {—il/) — m} . ( 10 - 2 ) 

®Note that for 0 < n < d/2 — 1 the difference — d/2 -|- n) — ipi^) is analytic at s = 0. 

^^Integration of Grassmann variables ip, ip satisfies J dip dip = a. 


19 





As a consequence, 

r = S' — logDet [ilp — m] = S — \ logDet 

= S+ \ r — e-™' Tr e-" . (10.3) 

2 Ja-2 t 

The heat-trace in the case of an electron of charge —e in the presence of a constant 
magnetic field B > 0 (Landau problem) in Euclidean four-dimensional spacetime 
can be readily computed from the spectrum of the squared Dirac operator ij) (see 
appendix [B]) , which is given by 

-|- 2ei? (n-|-1/2) ± ei? with A G n = 0,1,2,... (10-4) 


The two-dimensional momentum k corresponds to the timelike coordinate and the 
coordinate along the direction of the magnetic field. The second term in the eigen¬ 
values corresponds to a one-dimensional harmonic oscillator of frequency 2eB in 
one of the transverse directions, whereas the last term is the coupling between the 
magnetic moment of the electron and the external field. Apart from the degeneracy 
eil/27r per unit transversal area of these Landau levels there is an extra two-fold 
degeneracy corresponding to the electron/positron pair. From the spectrum (110.4p . 
taking into account the degeneracies, the heat-trace results 


Tr e 


_ 


eB 


= Vol(M^) — 


(fk 


—rk'^ 


2vr Jr 2 (27r)2 
,4 eB coth (eBr) 


OO 

2 ^ |g-'r2eB(n-|-l) g-T2eSn| 


n=0 


f ^4 eB 

jRi ^ 47r2 


(10.5) 


We therefore obtain an explicit expression for the one-loop corrections to the effective 
action 

r = ^[i?] -|- —j [ d^x [ e~™ coth {eBr) , (10.6) 

where the classical (Euclidean) action for the magnetic held reads 

5[.B] = [ d^x]-B‘^ . (10.7) 

jvk 2 


The divergence of the second term in expression (I10.6|) as A —)• oo can be removed 
upon the subtraction of the leading terms of the heat-trace for small r. These extra 
terms correspond to the renormalization of the cosmological constant and the electric 
charge. Let us write only the resulting hnite effective action after the subtraction of 
these counterterms. 


"f"*!=n L r $ 

Since the expression between braces is 0{t^) we could remove the cutoff A. The 
Hrst counterterm gives a held-independent contribution to the effective action so, as 
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mentioned, it is interpreted as a renormalization of the cosmological constant. On 
the other hand, the second term gives 


eB f 4 dr eBr 


i log {K^/w?) + .. .^ , 

(10.9) 


where the dots denote terms which are finite as A —>• oo. This contribution represents 
a divergent correction to the classical action (|10.7D . so it can be removed by a 
renormalization of the magnetic field, 

■Bphys ~ 247^2 B . (10.10) 


In order to mantain gauge invariance after quantum corrections the product of field 
and charge renormalizations must cancel to this order of perturbation theory; this 
gives the charge renormalization 

ephys = e ^1 - logA^ + .. .^ . (10.11) 


In consequence, the QED one-loop /3-function reads 


/3(e) = AdAe 


127r^ 


( 10 . 12 ) 


As in the case of the A cj)^ theory considered in section [51 the /3-function is positive. 
This implies that the coupling constant increases with the cutoff and could eventually 
become infinite for a finite value of A (Landau pole); equivalently, if the coupling 
constant e is kept finite then the renormalized charge Cphys vanishes signaling a 
complete quantum screening of the electric charge (triviality problem). Although 
the /3-function in this theory can only be computed perturbatively, there exists 
strong numerical evidence of the triviality of QED from nonperturbative lattice 
simulations. As a matter of fact, before 1970 the Landau pole problem seemed an 
unavoidable feature of any physically relevant field theory, until it was discovered 
that Yang-Mills theories have a negative /3-function [iniiiiiiiaiia]. We will carry 
out this calculation in section m 

Let us turn back to the finite expression (110.811 for the one-loop corrections to 
the effective action of the magnetic field. In Minkowski space, the effective action is 
related to the vacuum persistence amplitude 


e*'" = ( 0 ; out| 0 ; in) , 


so that 


|( 0 ;out| 0 ;in )|2 = . 


(10.13) 


(10.14) 


This implies that if there is an imaginary part in the effective action in Minkowski 
spacetime then the vacuum is unstable and the probability of pair creation is not 
zero. However, the magnetic field is unchanged under an inverse Wick rotation so 
the effective action (llO.Sp remains real in Minkowski spacetime and, as expected. 
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there is no particle creation in a purely magnetic field background. Nevertheless, the 
situation changes in the presence of a constant electric field. The full expression in 
the presence of both electric and magnetic constant fields is called Heisenberg-Euler 
effective action M- Since, in general, this action must depend on the invariants 
E‘2 _|_ ^2 E ■ B, we can obtain the effective action for an electron/positron 
system in the presence of solely a constant electric held by simply replacing B ^ E 
in expression (jlO.Sp . 

Now, upon an inverse Wick rotation, the electric held changes as E ^ —iE so 
(together with a global sign change) the effective action in a purely electric back¬ 
ground in Minkowski spacetime reads 

Although no imaginary contribution is manifest in this expression, it is certainly ill- 
dehned due to the singularities of the cotangent. Some prescription must be imposed 
in order to avoid these singularities; in other words, the integration path must be 
moved in the complex plane and thus imaginary contributions will appear leading to 
a nonvanishing probability of pair creation from the vacuunJ^. Since Euclidean time 
is dehned in the negative imaginary axis. Euclidean electric held corresponds to an 
analytic extension to the positive imaginary axis. Therefore, to recover the effective 
action in Minkowski spacetime we must approach the positive real axis from above, 
which means that the electric held in ()in.l5p has a small positive imaginary partly. 


r[E] 



|cot [e£'(l + zO)r] 


1 eEr 

J 



, .It 

cot (r -I- zO)-h - 

r o 


(10.16) 


This small imaginary part in the argument of the cotangent leads to a nonvanishing 
imaginary contribution due to the inhnitely many simple poles at r = vrn (with 
n = 1, 2 ,3,...) close to which 


1 


sin (r + zO) 


= I 


(-ir 


r — irn + zO 


= —(— l)”7r 5{t — zrn). 


(10.17) 


In consequence, the imaginary part of the effective action results 


I 



d^x 



(10.18) 


^^This imaginary contribution can also be understood from the fact that the series expansion in 
powers of eE is not Borel summable, as opposed to the purely magnetic case of expression (110.811 . 
which leads to an alternating series in powers of eB (see appendix m. 

^^Equivalently, since the change E —>■ —iE implies that eigenvalues that grow in the -|-oo-direction 
turn into eigenvalues growing in the —zoo-direction, one should introduce a positive imaginary part 
in the proper time r in order to have a well-defined heat-trace in the definition of the functional 
determinant. 
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This gives the pair production rate from vacuum induced by an external electric field, 
also known as Schwinger effect m- Expression (|10.18l) shows that the probability 
of pair creation is exponentially suppressed for electric fields below the critical value 

(mcj_ 1 ^13_lQi8^oit/m, (10.19) 

e he 

which is about three orders of magnitude above nowadays projected high-power laser 
facilities. Nevertheless, alternative theoretical settings which could enhance the rate 
of pair production are currently under study in order to estimate the possibility of 
reaching some experimental confirmation of this phenomenon in the next few years. 

11 /^-function in pure Yang-Mills 

In this section we show how by direct application of the heat-trace expansion 
(17.5p one obtains the /3-function in Yang-Mills theory. The non-abelian gauge field 
Afj_{x) = A‘^{x) Y“ takes values in the Lie algebra of some gauge group G with basis 
vectors (in some convenient representation) Y“, a = 1, 2,... , dimG. The invariant 
classical action in four-dimensional pure Yang-Mills is given by 

i / ‘r , (n.i) 

where g is the coupling constant. 


Ff,^ = i[D^,D^] 

(11.2) 

is the held strength F^^{x) = F“j^(x) Y“, and 


Df, = df,- iA^ . 

(11.3) 


In order to obtain the operator of quantum fluctuations, instead of computing the 
second functional derivative we make the shift A^[x) —)• A^{x) + a^(x) where 
now A^{x) is interpreted as a classical background field and the quantum gauge field 
is represented by a^{x) (Background Field Method). The second order expansion in 
af^{x) of the classical action then gives 

^ \ I tT{a^{-[Dty,[D„,af,]] - [Df„af,][D„,a„] + 2i , 

9 

(11.4) 


where the covariant derivative and the field strength are now computed 
exclusively from the background A^^. The functional integral that gives the effective 
action contains a divergence -due to the integration along gauge equivalent field 
configurations- which can be appropriately factorized by imposing a gauge condition 
on the quantum fields. If we choose the gauge condition = 0 and introduce 

in the quadratic part of the action the gauge fixing term in the Feynman gauge 
(^ = 1), then the second term in the r.h.s. of (|11.4I) is cancelled and we get 


c(2) 

gauge 



(11.5) 
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where the operator of quantum fluctuations of the gauge field is now given by 

-J's-gauge = - [Dp. •]] + [F^u, •] = -V (5 - F^^ . (11.6) 

On the other hand, after imposing the gauge choice = 0, gauge invariance 

is ensured by the introduction of (Grassmann) ghost fields with action 

'S'ghost = / dx tr (-c(x)[i:)^, [Il^,c(x)]]) . (11.7) 

The corresponding operator of quantum fluctuations of the ghost held then reads 

d'S-ghost = -[Dp, [Dp, •]] = -(a - . (11.8) 


Now, the one-loop effective action can be expressed in terms of the functional de¬ 
terminant of these operators, 

r[^] = S'[A] -I- ^ logDet {(5^5'gauge} - logDet [ 5 ^ *S*ghost } • (11-9) 

There is no — ^ factor in the third term (as opposed to the second one) because, as 
already mentioned, the auxiliary ghosts are Grassmann helds (see also section □TO]). 
The functional determinants in expression (Ill.Qp can be regularized by means of the 
heat-traces of the corresponding quantum fluctuation operators. 


1 1*00 J . , 

r[^] = S[A] - - / — ( Tr e-^'^'^eauge _ 2Tr _ (n ig) 

2 7^-2 r V / 


Notice that we have also introduced an IR regulator m. If we replace in this expres¬ 
sion the corresponding heat-trace expansions, the only divergent contributions arise 
from the terms proportional to the 02 coefficients. According to (|7.9h . for the gauge 
held, 

^gauge ^ 1 4C2 (g) ^ 4C2(g) ^ C2{g) , (ii.ii) 

whereas for the ghost held, 

^ghost ^ (^^^^2) 


The Casimir factor C 2 {G) arises from the trace trA"“A"^ = G 2 {G) 6 ab in the adjoint 
representation. Replacing these contributions to the heat-traces into (I11.10|) we 
obtain 


r[A] = 5[A] 

= 5[A] 


- f —e [ dx — 2 02 ^“*)-I-(hnite terms) 

2 7a-2 t (dvr)"' 7 ]r4 2 

C' 2 (G) log (A/m) [ dx -h (hnite terms) . (11.13) 


As is well-known, UV-divergencies in gauge theory are logarithmic. The divergence 
in (I11.13|) can be removed by a redehnition of the coupling constant g that appears 
in the classical action S'fA], 


5phys 
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247r2 


G 2 {G)g‘^ log(A/m) 


(11.14) 
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The /3-function then reads 


/3 = Kd^g = C2{G) . (11.15) 

Note that, being /3 < 0, if the coupling constant is small at large energies then it 
decreases as A —>■ oo so that there is no Landau pole and perturbative calculations 
become appropriate for describing short-distance processes. This behaviour of the 
coupling constant is consistent with the anti-screening effect of virtual gluons. 


12 Quantum Field Theory at finite temperature 


The partition function Z[J] given by expression (18.ip can also be used to determine 
the thermodynamical properties of quantum fields at finite temperature. The stan¬ 
dard procedure consists in performing an analytic continuation to Euclidean time t 
in the operator of quantum fluctuations and imposing periodic boundary conditions 
at t = 0 and t = (5, where /3 stands for the inverse temperature. In order to illustrate 
this method we consider the simplest example, namely, a non-interacting scalar par¬ 
ticle of mass m conhned in a one-dimensional box of length L. After implementing 
the Wick rotation, the partition function (for J(x) = 0) reads 

log Z = — i log Det {—(9^^ — (9^ -|- I 

= ^C'(0) + ^C(0)log/r2, (12.1) 

where the ^-function is given by 



1 

f(^ 





EE' 

n=l /cGZ 


4-77^ 


( 12 . 2 ) 


Note that the first term in the eigenvalues corresponds to the oscillating modes 
along the compactified Euclidean timelike direction (Matsubara frequencies). If we 
consider fermionic particles, we should impose antiperiodic boundary conditions at 
t = 0 and t = (5, instead. The second term in the eigenvalues comes from imposing 
Dirichlet boundary conditions at a; = 0 and x = L. In the second line of (112. 2p 
we have applied, as usual, the Mellin transform that relates the (^-function to the 
heat-trace. 

Next, we use Poisson inversion formula in the sum over /c G Z and afterwards, 
for the resulting terms corresponding to A: = 0, we use Poisson inversion in the sum 
over n = 1,2,3,... After all these manipulations, the result reads 


C(s) = - 




4i/7rr(s) 

/3L r 
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47rr(s^ 

°° 2 

n=l 


f 
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(12.3) 
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Finally, we perform all integrations in r to get 


^ /3mi-2"r(s-i) ^ /3Lm2-2" ^ ^ ii:i_,(2Lmn) 

“ A /z:-nf„\ ^ /!_ ! „ TT rF73“ ^ 

n=l 


+ 22 


4y^r(s 

i-s 


47r (s — 1) ' 7rr(s) * 


-i/ttFIs) 

'' ' n=l fc=l 




TT^re^ ^ 4 2 


L2 


If 2 ii'i_ I3k\ rrf + 


'K’^rf 

~1X 


(12.4) 


Due to the 1/F(s) factor, the function C('S) at s = 0 is given exclusively by the 
second term 

C(0) = 


47r 


(12.5) 


On the other hand, the derivative C{s) at s = 0 readJ^ 


, , mB jBLm? , , ns Bm Ki(2Lmn) 

C O = ^ + ^1-— -i + iogm^) + ^ y -1 + 

^ ^ 2 dvr ^ ^ vr ^ n 

n=l 


+2EE 

n=l k=l 


.-Pk^r, 


( 12 . 6 ) 


We can now use expression (|12.1I) to evaluate the partition functiorP^. 

‘“gz = ^ (-1 + log(™V m^)) + ^ E^4^ 

-Eiog(i-oW 


+ 


n=l 


n=l 


L2 


(12.7) 


This expression depends on the arbitrary scale /r -due to the nonvanishing value of 
C(0)“ and must be thus appropriately redefined by means of some renormalization 
prescription, as the vanishing of logZ for inhnite values of m?. This determines 
the partition function for a massive scalar field confined in the interval [0,L]. For 
simplicity, we consider next the massless case 




n=l 


— e 


3 Trn. 
^ L 


Note that the first term is proportional to the vacuum energy 


n=l 




“ ~ ~1al 


( 12 . 8 ) 


(12.9) 


s=-l 


whereas the second term is the partition function for the grand canonical ensemble. 
We can now compute, at temperature /3, the mean energy 


E-= + 
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24L ' L ^ - 1 ’ 
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(12.10) 




= -log(l-ir). 
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Figure 2: Mean energy, pressure and entropy for a scalar massless particle in a 
one-dimensional box of unit length. 

the pressure 

- OO 

p := - log Z = + J 2 J 2 e/3^ _ 1 ’ 

n=l 

and the entropy 

OO OO 

S’ := log Z -h /3 F; = - ^ log (^1 - Z ^ —l ' (12-12) 

n=l n=l 

These thermodynamic quantities are represented in figure [2] as a function of the 
inverse temperature (5. At high temperatures (small (3) all these three quantities 
diverge, whereas at zero temperature (/3 —)• 00 ) the entropy vanishes, as expected, 
but the pressure and the mean energy attain a finite non-zero value corresponding 
to the vacuum oscillations of the quantum field. 

13 Noncommutative quantum fields 

Our last application of functional methods in quantum field theory is related to the 
computation of n-point functions from the Seeley-de Witt coefficients. In particular, 
we will compute one-loop corrections to the 2-point function in a noncommutative 
field theory in order to show the so-called UV/IR mixing effect. 

Noncommutative field theories are formulated on a spacetime whose coordinates 
satisfy 


— 2z , (13.1) 

where 0 is an antisymmetric matrix with dimensions of squared length. The 
fields are thus operators </>(x) that depend on these noncommutative coordinates. 
By means of the Weyl-Wigner transformation this space of noncommutative fields 
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is isomorphic to the algebra of ordinary functions (/>(x) (defined on the usual com¬ 
mutative M'^) with the following noncommutative ^-product: 

ct){x)-kx{x) (13.2) 


Under this nonlocal associative product, know as Moyal-product, ordinary coordi¬ 
nates satisfy Xi-kXj — Xj-kXi = 2i Qij. We will consider the case where 0 has maximal 
rank and that, after an appropriate choice of coordinates, can be written as 
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-9 

0 / 


(13.3) 


9 G M'*' is now the parameter that characterizes noncommutativity. In this formu¬ 
lation the Euclidean action of a noncommutative scalar field with a quartic self¬ 
interaction reads 

S[4>]= [ dx (1)^ + ^ , (13.4) 

where (j)^ := (p-kcp-kcp-kfj). Since the ^-product and the ordinary commutative product 
differ in total derivatives, quadratic terms in the action do not need to be modified 
in Moyal space. 

After computing the second functional derivative, the field fluctuation operator 
results 

6^S = -52 + ^ {L((/.2) + + L{cj))R{cP)} , (13.5) 

where L and R represent (left- and right-) Moyal-multiplication. In consequence, 
the one-loop effective action is given by 

r[</,] = 5[0] _ 1 ^ Tr e'" (13 5 ) 

2 Jo X 

Since the field fluctuation operator is nonlocal, the standard heat-trace asymptotics 
described in section [7] no longer holds HZj. Let us study the effects of nonlocality 
on the quantum propagator (or 2-point function), which is given by the part of the 
effective action that is quadratic in the field cj). 

As can be seen from (113.6p . there are three contributions to the effective action 
which are quadratic in the field arising from each term L{(j)‘^), R{(t)^) or L{(j))R{(j)) 
in an 0(A) expansion of the exponential. The contribution of the first two of these 
terms is the same as in the commutative case: it is a divergent quantity that is 
regularized upon a mass renormalization. The most remarkable consequence of 
nonlocality on the heat-trace is due to the last term in the operator 6‘^S, which 
mixes left- and right- Moyal multiplication. As we will see, the contribution of the 
term L{4>)R{(j)), though UV-finite, leads to an IR-divergence in higher-order quantum 
corrections. In the Feynman diagram formulation of noncommutative field theories 
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this term corresponds to a one-vertex nonplanar diagram. In order to evaluate this 
contribution we need the following result: 


Tr 
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(dvrr)^ 
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/ dx 4 >{x) e 
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(f){x) 
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--rP 


(13.7) 


(4vrr)2 7 ]r 4 (27r)4 

which has been obtained in m using the worldline formalism. In the second line we 
have used 0^ = —0^ I 4 . In consequence, the nonplanar contribution to the quadratic 
effective action reads 

.,2 




iIo.m 2 Ki{2me\p\) 

\4 l^(p)l 


(27r)^dSvr^ 

This expression shows that the one-loop propagator contains a contribution 


(13.8) 


(13.9) 


which is finite but diverges for small values of the momentum, even for massive helds. 
As a consequence, although this contribution is UV-finite, the use of this propagator 
in higher-order calculations (where p represents the momentum of virtual particles 
and is thus off-shell) introduces IR-divergencies. This phenomenon, known as UV/IR 
mixing [IS], can be interpreted in the following way. Note that the ^-product can 
be written as 


(<j)-k <j))(x) = -—— [ dx'dx" (f)(x')(f)(x 
det 0 Jks 


^-10^2 {x'-x)i(x"-x)j 


(13.10) 


Therefore, if the support of (p{x) is restricted to a distance di in the xi-direction, 
then (j)l is nonzero in a distance 62 ~ O /61 in the a: 2 -direction. This means that 
a wave packet of width di <C Vd gives a function (pl of width 62 ^ Vd- In other 
words, small pulses are largely spread upon self-interaction, even for massive fields 
-which have exponentially decreasing propagators. As a consequence, scattering at 
low energies is affected by high energy virtual particles. 
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A Divergent sums 

In this appendix we will illustrate with examples some criteria which are used to 
define divergent sums. In order to make contact with analytic extensions it is prefer¬ 
able to work with alternating sums, so we define 

r7R(s) = 1-* - 2-* + 3 -* - 4 -* + ... (A.l) 
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if lZ{s) > 1, or as its analytic extension otherwise. Its relation with the Riemann 
C-fnnction can be derived from the expression 


1_, _ 2-* + 3-* _ ,,, = ( 1 -^ + 2-" + 3-" + ...)- 

-2 X + + + (A.2) 

which, extended to the whole complex plane, gives 

i1r{s) = (1 - 2^“^) C,r{s) . (A.3) 

From the values of Cij('S) at s = 0,-1,—2 -computed in section [3}- we get the 
following values of the 77 -function: 

77^j(0) = i = l- l + l- ... (A.4) 

= i = 1-2 + 3-... (A.5) 

77^j(-2) = 0 = 1-4 + 9-... (A.6) 


We will next compare this definition of the alternating divergent sums in terms of 
the analytic extension of the r 7 -function with some well-known criteria. 

Let us first consider Cesaro summation which, if the sequence of partial sums 
does not converge, defines the divergent series as the limit, if it exists, of the averages 
of the divergent sums. Take, for example, the series 1 — 1 + 1 — ... whose partial 
sums give 1,0,1,0,... Although this oscillating series is not convergent the averages 
of the partial sums give the sequence that converges to Therefore, 

we write 

l-l + l-... = i (C,l), (A.7) 

where (C*, 1) stands for “the first average of partial sums”. To explain this notation, 
let us next consider the series 1 — 2 + 3 — ... whose partial sums give 1,—1,2,—2,.... 
The sequence of averages of this sequence is 1,0, |, 0, |, 0,... which oscillates around 
the two values 0 and However, a new average of these averages gives the sequence 
1, . which converges to Therefore, we write 

1-2 + 3-... = i (C, 2 ), (A. 8 ) 

where [C, 2 ) stands for the second average which was required to consider in order 
to obtain a convergent sequence. Following the same arguments one can prove that 

l-4 + 9-... = 0 (C',3). (A.9) 

As expected, the more divergent the absolute values of the term in the series are, 
the more number of averages have to be taken into account. The interpretation of 
this method is that if a series X] is not summable, then the sequence ai,ai + 
02 , oi + 02 + 03 ,... is not convergent. The average of partial sums instead give the 
sequence oi, oi + ^ 02 , oi + | 02 + | 03 ,... The idea behind (C, 1) is that instead of 
adding a full term a-n at each step of a partial sum, only a fraction of it is added; 
of course, this fraction becomes closer to one as the number of steps increases. As 
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one can easily verify, a higher number of averages means that the fraction of a given 
term which is added at each partial sum approaches 1 slower. 

It can be proved that (C, n + 1) is stronger than (C, n), meaning that if a series 
can be summed according to (C, n) then it can also be summed (to the same value) 
according to (C, n + 1). In general, Cesaro summation is linear in the series, regular 
(it is stronger than the usual criterion for convergent series) and stable (if one 
removes the first term one gets a summable series which differs from the original 
one in this term). 

Let us next consider Abel criterion, which defines the value of a series ^ an as 

CXD C» 

an'■= lim Onx", (A.10) 

n=\ n=l 


as long as the r.h.s. converges for x < 1 and the limit exists. Note that this method 
coincides with a heat-kernel type regularization -based on the cutoff function if 

we define x := e“L It can be shown that Abel criterion is linear, stable and stronger 
than Cesaro criterion. For example, taking the limit x —)• 1“ in the convergent 
expansions 

1 — X + x^ — ... (A.11) 

1 - 2x + Sx^ - ... (A.12) 

1 - 4x + Ox^ - ... (A.13) 

one obtains the same values as in ()A.7p . (lA.Sp and (jA.Op . 

Finally, let us consider Borel summation, which is stronger than Abel criterion. 
To understand Borel summation it is convenient to write 


1 

1 + X 

1 

(1 -|- x )2 

1 — X 
(1 + x)3 


CX) 

n=l 



(A.14) 


and then formally interchange the sum with the integral to define 



CX) 



n=l 


(A.15) 


One can verify that Borel summation also gives the values in ()A.7p , (IA. 8 D and (|A.9I) . 
The strength of this method can be understood in the divergent series X^(—I)"" nl 5 "', 
for positive g. Borel summation readily gives 


^(-l)-n!ff- 

n=l 


1 , e-^/9 

- dr - 

g Jo 1 + T 


(A.I 6 ) 


which leads to the result 0! —1!-|-2! —... = 0.596347... The importance of expression 
(IA.16P is that it indicates that a divergent series which is Borel summable might 
occur in a physical problem as the result of expanding a non-perturbative expression. 
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In general, a series is Borel summable only if ^anT^/n! converges 

in some neighbourhood of the origin to some function B{t) that can be analyti¬ 
cally extended to the positive real axis. Then is defined as the integral 

g-^Jo^e-/9B{T)dT, as long as this integral exists. In fact, under certain assump¬ 
tions on analyticity and on the decrease of the remainder, Watson’s theorem states 
that an asymptotic series uniquely determines the sum in terms of an integral of 
Borel’s type m- 

It it interesting to apply this definition to expressions (110.811 and (I10.15p for the 
effective action of an electron/positron system in constant magnetic and electric 
backgrounds, respectively. In the purely magnetic case, the function B{t) is given 
by (cothr — l/r — r/3)/r^, which is analytic in a neighbourhood of the positive 
real axis and -upon integration in r- leads to an alternating series (in powers of 
g = eB/m?) for the effective action which is Borel summable. On the contrary, for 
the purely electric case, B{t) is given by the function (cot t — Ijr -\- t , which 
has infinitely many simple poles in the positive real axis and, if integrated in r, leads 
instead to a divergent series (in powers oi g = eE/m?) of positive terms. In the 
presence of singularities some prescription must be implemented, as was shown in 
section uni in order to perform the r-integral [20] . 


B Landau problem 


Let us consider a massless Dirac spinor of charge — e in the presence of a constant 
magnetic field B in the positive z-direction. The corresponding Dirac operator is 


-0 = - ieAfj) , 


(B.l) 


with = (0,0,Ba:,0). In our notation = {t,x,y,z). We use the metric = 

(H-) and, as a representation of the Clifford algebra { 7 ^, 7 ^} = 2 rj^'^, the Weyl 

basis: 
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0 

I2 




The squared operator 


= {d- ieAf - 5 ( 7 ^, 7 ''] ieF^^ 
= □ — 2ieBxdy -(- e^B^x^ + eB 


0 \ 

0 ) 


restricted to the subspace 

^{t,x,y,z) 


p{x) 

X{x) 


^—ikot—ik2y—ikzz 


(B.2) 


(B.3) 


(B.4) 


can be written as 

U)^ = -kl + kl-dl + e^B^ix - k^jeBf + ^3 ) • (B.5) 
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Therefore, in Euclidean spacetime {ko —>■ iko) the spectrum is 


fcg + + 2eB{n + 1/2) ± eB 


(B.6) 


with ko, ^3 G M, n = 0,1, 2 ,... and an extra double degeneracy corresponding to the 
two Weyl spinors. There is also an infinite degeneracy because the eigenvalues do not 
depend on k 2 - However, we can measure this degeneracy in terms of the area LxLy 
transversal to the magnetic field. The shift k 2 /eB in the x-direction imposes the 
constraint \k 2 \/eB < Lxl2. Since the number of states corresponding to this range 
of momenta is 2 x eB x Lxl2 x Lyl2'K we obtain a total degeneracy LxLy x eB/2 'k. 

Let us next consider the massless Dirac spinor of charge —e in the presence of 
a constant electric field in the x-direction. Using = (—Ex,0,0,0) for the Dirac 
operator defined in (IB.lh we obtain 


.0^ = 0 + 2ieExdo — e^E^x^ + ieE 


0 

0 - 0-1 


which restricted to the subspace 
T(t,x,y,^;) 


4>{x) 

X{x) 


^—ikQt—ik2y—iksz 


reads 


Ij)'^ = kl + kl- dl- e^E‘^{x - ko/eEf + ieE 


0-1 0 \ 

0 -uV • 


(B.7) 


(B.8) 


(B.9) 


Note that the operator contains an inverted harmonic oscillator. Tunneling across 
this barrier gives the pair production rate. Alternatively, the introduction of a 
positive imaginary part in the proper time provides the exponential decrease that 
allows the use of heat-kernel techniques. 
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